Second-Order Partial Differential Equations 

Lecture 4 

Second-Order Partial Differential Equations 



In this lecture and in the next, we'll briefly review second-order PDEs. We'll begin with one of 
the simplest of such PDEs: the Laplace equation. 

A The Laplace Equation 

The Laplace equation appears in many branches of physical sciences, two of which being 
electrostatics and fluid mechanics. 

Let the Laplace equation hold in the upper-half plane in a two dimensional space: 

= ^^y > 0- (4-1) 



Needless to say, there are (infinitely) many solutions to this equation, corresponding to (infinitely) 
many different physical situations in the realms of static electricity or incompressible and irrotational 
fluid flow. In order for the solution of this equation to be unique, we need to impose some conditions 
on u. 

Consider the square in the upper-half plane with the x-axis between -R to R being its lower 
horizontal boundary. We will think of the upper-half plane as this square in the limit R ->■ co. We 
require that the value of u is equal to zero at all sides of the square except the lower horizontal side. 
In the limit R ^ x>, this means that 

^mu{x,y) = 0, (4.2) 

where 



In addition, we require that, on the the lower horizontal side of the square which is the x-axis, 

u{x,0)=Ax). (4.3) 

Equations (4.2) and (4.3) constitute a set of boundary conditions for u. More precisely, the values of 
u on three sides of the square are required to be zero, while those on the fourth side are given hyf{x). 
Equation (4.1) together with these boundary conditions for m is a boundary-value problem. 

It may help our perspective if we compare this with the much simpler problem of the Laplace 
equation in the one-dimensional space, which is 

^uix) = 0. 

Being of the second order, the equation gives the solution 

m(x) ^ ax + h, 

which has two arbitrary constants. To determine the constants a and b, two conditions on u are 
required. 

We may choose these two conditions to be at the same point of x. For example, we may choose 
the two conditions to be the value of u and the value of m' at x = 0, and the problem is to find the 
solution for x > 0. Such a problem is called an initial-value problem, and the solution is uniquely 
given by M = m(0) -I- u{0)x. 

Alternately, let the range of x for the equation be [0, 1]. Then we may choose these two 
conditions to be the values of u at the endpoints, say m(0) = 1 and m(1) = 2, which are conditions 
imposed at two different points.This is called a boundary-value problem. The solution of this 
problem is uniquely given by w = x -i- 1. 

Let us see what happens if, in a boundary- value problem, we specify the values of u' instead of 
those of u. We see immediately that there is a difficulty. For, if m is a solution of u" = 0, u' is a 
constant. Therefore, there is no solution satisfying these boundary conditions unless the two values 
of m' specified at the endpoints are equal. But there is another difficulty if the specified values of u' 
at the two endpoints are equal. This is because if m is a solution of this boundary-value problem, so is 
u + c, where c is an arbitrary constant (why?). Thus the solution of such a boundary- value problem is 

— 1 — 



Second-Order Partial Differential Equations 
not unique. 

There are also other kinds of boundary conditions. For example, we may specify the value of u at 
one of the boundary points, and the value of u' at the other boundary point. As an example, let 
m'(0) = 1 and m(1) = 0. Then the solution of the boundary-value problem exists and is unique, and 
is given by u{x) = (x- 1). 

Let us return to the two-dimensional Laplace equation (4.1) with the boundary conditions (4.2) 
and (4.3), which specify the values of u at the boundary of the upper-half plane. To obtain the 
solution to this boundary- value problem, we focus on the x-dependence of u{x,y) on x and denote the 
Fourier integral of u in the x variable by 



u(k,y) = J 



e-'^u{x,y)dx. (4.4) 



This Fourier transform exists because u(±co,y) = 0. Instead of solving for u{x,y), we shall solve for 
u{k,y). Once u{k,y) is found, we immediately obtain u{x,y) with the use of the Fourier inversion 
formula. 

Problem for the Reader: 
Substitute (4.5) into (4.1), obtain an ODE for u{k,y) and solve this equation. Use the boundary 
conditions to determine u{k,y). 

Answer 

Since D operating on e''"' is equal to ik, where D = we get 



-k^ + ]uik,y) = 0. (4.6) 



,0 

u{x,y) = J 



dy^ 

We see that the PDE becomes an ODE, which is easily solved. This is the advantage of expressing u 
as a Fourier integral. Solving the ODE, we have 

uik,y) =Aik)e''y + B(k)e^'^. (4.7) 

The boundary value condition (4.2) forbids the solution to blow up as y oo. Thus we may eliminate 
one of the constants in the solution (4.7). The condition (4.3) enables us to find that 

^e"yf{k), (4.12) 
In 

which satisfies both the Laplace equation (4.1) and the boundary conditions (4.2) and (4.3). In the 
above, /(^) is the Fourier transform off{x). 

The factor e""^'^ provides a convergent factor for the Fourier integral in (4.12) when y is non-zero. 
Thus the integral in (4.12) converges nicely as one leaves the x-axis, where the values of u are given, 
into the x-y plane, where y non-zero. The boundary- value problem of the Laplace equation is 
therefore called well-posed. It is also the kind of problem that people naturally ask themselves, an 
example being "what is the static electric potential inside Benjamin Franklin's silver pint can?" 

As we recall, for the first-order PDE (2.15), all values of m on a characteristic curve are 
determined once an initial value on this curve is given. For the Laplace equation, one needs the 
values of u at all points of the boundary to determine the value u{x,y) at any point in the interior 
region. Trivially, no value of the solution at any interior point is known if the value of u is given only 
at one point. 

Problem for the Reader: 

If/(^) = , ^ 2 > find u{x,y). 
Answer 



"ated out ; 
u{x,y) = 



■lifei 



The integral above can be easily integrated out and we get 

x^ + (1 +y)'^ ' 

We next consider the initial- value problem of the Laplace equation. More precisely, let u 
satisfying the initial conditions 

m(x,0) =j{x), 

and 
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Uy{X,0) = gix). 

We also require that 

M(±oo,y) = 0. 

The last condition ensures that u(x,y), considered as a function of x, has a Fourier transform. 

As before, we think of the upper-half plane as an infinite square in the upper-half plane. In the 
problem we shall treat, two conditions are given on the lower horizontal boundary of the square, one 
condition is given to each of the two vertical boundary lines, but no condition is imposed on the 
upper horizontal boundary line. Therefore, this is, more precisely, an initial- value problem with 
respect to the y variable and a boundary- value problem with respect to the x variable. 
Problem for the Reader: 
Solve the Laplace equation (4.1) with the initial conditions given above. 
Answer 

sinh ky 



fik) cosh^j + gik)- 



(4.17) 



k 

Since both cosh ky and sinh ky blow up exponentially as a function of /c as goes to infinity, it is 
impossible for the integral in (4.17) to converge for all values of y unless both /(A:) and g{k) vanish 
more rapidly than an exponential function as \k\ goes to infinity. 

We call the initial-value problem of the Laplace equation not well-posed. 

B The Wave Equation 

Consider the wave equation 

^xx-^tt = 0. (4.24) 

The wave equation describes many physical phenomena, a notable example being the propagation of 
the electromagnetic wave in the vacuum. 
Problem for the Reader: 

Let 

u = x + t,v = x-t. (4.25) 
Express the wave equation (4.24) with the independent variables u and v. 
Answer 

dudv 

From the equation above, we see that the solution of the wave equation is a sum of an arbitrary 
function of u and an arbitrary function of v: 

^(x, t) = f{u) + giv) = fix + t) + g{x - t). (4.26) 

We note from (4.26) that the solution of the second-order wave equation has two arbitrary functions. 
This is to be compared with that of a first-order PDE of two independent variables, which generally 
has one arbitrary function. We also note ihaXfix -i- f) is a wave travelling to the left with unit velocity, 
and g{x - ?) is a wave travelling to the right with unit velocity. 
Problem for the Reader: 
Let ^ satisfy the following initial conditions on the x-axis: 

^(x,0) = a{x),^t{x,0) = b{x). (4.27) 

Find (l>{x,t). 
Answer 

It is straightforward to show from (4.26) and (4.27) that 

. a{x + t) + a{x-t) 1 f-^^S ^ i\ 1 I /A 
(j>{x,t) = ^2"-^ ^ + ^ J b{x')dx'. (4.29) 

Obviously, we may also impose the Cauchy data on other lines or curves, not just the x-axis. 

See the textbook for more detailed explanations, more examples, as well as the treatment of 
other partial differential equations. 
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